Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN WF$t style file v2.2) 



m 

O 
O 

(N 



MS 1603.6+2600: an Accretion Disc Corona source? 

P.G. Jonker 1 *, M. van der Klis 2 , C. Kouveliotou 3 , M. Mendez 4 , W.H.G. Lewin 5 , 
T. Belloni 6 

1 Institute of Astronomy, Madingley Road, CB3 OH A, Cambridge 

2 Astronomical Institute "Anton Pannekoek", University of Amsterdam, Kruislaan 403, 1098 SJ Amsterdam 

3 SD 50, Space Science Research Center, National Space Science and Technology Center, 320, Sparkman Drive, Huntsville, AL 35805, 
Nasa's Marshall Space Flight Center 

4 SRON, National Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands 

5 Department of Physics and Center for Space Research, Massachusetts Institute of Technology, Cambridge, MA 02138 

6 INAF, Osservatorio Astronomico di Brera, Via E. Bianchi 46, 23807 Merate, Italy 

2 February 2008 

ABSTRACT 

We have observed the eclipsing low-mass X-ray binary MS 1603.6+2600 with Chandra 
for 7 ksec. The X-ray spectrum is well fit with a single absorbed powerlaw with 
an index of ~2. We find a clear sinusoidal modulation in the X-ray lightcurve with 
a period of 1.7±0.2 hours, consistent with the period of 1.85 hours found before. 
However, no (partial) eclipses were found. We argue that if the X-ray flare observed 
in earlier X-ray observations was a type I X-ray burst then the source can only 
be an Accretion Disc Corona source at a distance between ~ 11-24 kpc (implying 
a height above the Galactic disk of ~8-17 kpc). It has also been proposed in the 
literature that MS 1603.76+2600 is a dipper at ~75 kpc. We argue that in this dipper 
scenario the observed optical properties of MS 1603.6+2600 are difficult to reconcile 
with the optical properties one would expect on the basis of comparisons with other 
high inclination low-mass X-ray binaries, unless the X-ray flare was not a type I 
X-ray burst. In that case the source can be a nearby soft X-ray transient accreting at 
a quiescent rate as was proposed by Hakala et al. (1998) or a high inclination source 
at ~15-20 kpc. 
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1 INTRODUCTION 

Low-mass X-ray binaries (LMXBs) are binary systems in 
which a <; 1 M© star transfers matter to a neutron star or 
a black hole. A large fraction of these LMXBs are found 
in the Galactic Bulge (see |van Paradijs fc White 1995| l. 
MS 1603.6+2600 was discovered with the Einstein satellite 
as a faint source at a flux level of ~ 1 x 10 _12 ergcm -2 s^ 1 
(0.3-3.5 keV IGioia et al. 19901 . |Morris et al. (1990)| found 
an optical counterpart that shows partial eclipses with a pe- 
riod of 111 minutes; they also reported that the depth of 
the optical eclipse is anti-correlated with the optical lumi- 
nosity. No radio emission was detected with an upper limit 
of 0.3 mjy (IMorris et al. 19901 . RO "SAT and ASCA detected 
this source at flux levels of ~ 1 x 10~ 12 erg cm -2 s _1 (0.1- 
10 ke V) and ~ 4 x 10~ 12 erg cm" 2 s" 1 (0.7-10 keV), respec- 
tively JMorris et al. 19901 we determined the flux from the 
spectral properties and count rate measured with ROSAT; 
IHakala et al. 19981 . 



It was immediately realized that this system is either a 
Cataclysmic Variable (the compact object is a white dwarf) 
or an LMXB (the compact object is a neutron star or a black 
hole) and that if this system is an LMXB its distance may 
be large (IMorris et al. 199011 . However, |Hakala et al. (1998)| 
proposed that the system could be a nearby soft X-ray 
transient in quiescence accreting at a low rate. They also 
showed that the system properties do not fit any of the Cat- 
aclysmic Variable categories. [Ergma fc Vilhu (1993)| showed 
that there are three possible evolutionary scenarios for 
this system and they argue that the system contains a 
neutron star rather than a white dwarf compact object. 
|Mukai et al. (2001)| detected an X-ray flare, which could be 
a type I X-ray burst. If the flare was indeed a type I X-ray 
burst it establishes the nature of the compact object as a 
neutron star. 
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In this Paper we report on a Chandra observation of 
MS 1603.6+2600. 
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2 OBSERVATIONS AND ANALYSIS 



The Chandra satellite fWeisskopf 1988^ observed 
MS 1603.6+2600 with the Advanced CCD Imaging 
Spectrometer (ACIS) on March 11, 2002 for nearly 7 ksec 
starting at 08:19 (TT; MJD 52344.346). The data were 
processed by the Chandra X-ray Center; events with ASCA 
grades of 1, 5, 7, cosmic rays, hot pixels, and events close to 
CCD node boundaries were rejected. We used the standard 
CIAO software to reduce the data (version 2.3 and CALDB 
version 2.21). A very weak streak caused by the arrival 
of photons during the CCD readout period (40 ^seconds) 
was present. We removed this streak before processing the 
data further using the ciao tool ACISREADCORR taking into 
account that we had windowed the CCD to 1/4 to reduce 
effects of pile-up. 

We detect only one source. After applying the ciao 
web-based tool fix offsets to correct for known as- 
pect offsets, we derive the following coordinates for 
the source: R.A.=16h05m45.88s, Decl -+25°51'45.1" (typ- 
ical error .6", equinox 200 0.0). The USNO B1.0 
coordinates iMonet et al. 200*31 of the optical coun- 
terpart of MS 1603.6+2600 [UW Coronae Borealis; 
USNO B1.0 Id. 1158-0232558, R.A -16h05m45.868(3)s 
Decl -+25°51'45. 560(3)"] are fully consistent with this. 

The source is detected at a count rate of 0.72+0.01 
counts per second; this could yield a pile-up fraction of ap- 
proximately 25-30 per cent for our frame time of 0.84 s. In 
order to investigate this further, we study the radial pro- 
file of the point spread function of the source (see Figure [*f}. 
Comparing the profile (Figure[""J with the theoretical profiles 
shown in figure 6.24 of the Chandra Proposers' Observatory 
Guide v.5 1 we note that the number of counts per pixel is 
lower in the inner 1-2 pixels than it would have been in the 
absence of pile-up. 

The spectra are extracted with 20 counts per bin. We 
only include energies above 0.3 and below 8 keV in our spec- 
tral analysis since the ACIS timed exposure mode spectral 
response is not well calibrated below 0.3 keV and above 8 
keV. We fit the spectra using XSPEC HArnaud 1 996) ver- 
sion 11.2.0ao, including an extra multiplicative model com- 
ponent in the fit-function, the ACISABS model 2 , to correct 
for additional absorption due to contamination by the opti- 
cal blocking filters in all our spectral fits. Since this model is 
only accurate to approximately 10 per cent we furthermore 
include a 10 per cent systematic uncertainty to the channels 
below 1 keV (channels 1-69). To correctly model the effect 
of pile-up on our observed spectrum, we fitted the spectrum 
using the pile-up model of |Davis (2001)| 

The spectrum is not well fit by a single absorbed black- 
body component [reduced x 2 = 2.9 for 173 degrees of free- 
dom (d.o.f.)]. A single absorbed power law yields a good 
fit, with reduced \ 2 = 1.0 for 173 d.o.f. (see Figure 
The power law index is 2.0+0.1 for an interstellar absorp- 
tion, N H , of (1.5 + 0.2) x 10 21 cm~ 2 . This N H is incon sistent 
with the value derived by |Dickey fc Lockman (1990)|(^ 4 x 
10 20 cm~ 2 ). However, using ASCA data |Mukai et al. (2001)] 
find that the absorbing column to this source varies as 
a function of the binary orbital phase; it can be as 



1 av ailable at |http://asc. harvard.edu/udocs/docs/docs. html| 

2 see http://asc.harvard.edu/cal/Acis/Cal_prods/qeDeg/ 



high as 4.6xl0 21 cm" 2 . The unabsorbed flux (0.1-10 keV) 
is 1.2xl0 -11 ergcm _2 s _1 . We present the parameters 
for this model in Table ["fl To test the consistency of 
the XSPEC pile-up model implementation we also fit- 
ted the spectrum with the ISIS package version 1.1.3. 
( Houck & Dcnicola 2000 ) and their implementation of the 
pile-up model of |Davis (2001)| We find that the results are 
consistent within the 90 per cent uncertainties. We note that 
in the ISIS fits we corrected the auxiliary response file using 
the CIAO tool corrarf to take the degredation of the ACIS 
quantum efficiency due to the contamination of the optical 
blocking filters into account. We present the ISIS fit-results 
in the second line of Tabled 

We create lightcurves with a time resolution of 10 and 
200 seconds using photons with energies 0.3-10 keV ex- 
cluding a circular area of a 2 pixel radius centred on the 
source position to mitigate the effect the pile-up has on 
the lightcurve. The lightcurve with a resolution of 10 s 
was searched for the presence of type I X-ray bursts, but 
none was found. In Figure [3] we plot the lightcurve in the 
0.3-10 keV energy band. We could not phase connect our 
observation with previous observations since due to the 
uncertainties in the ephemeris, and the long period be- 
tween our observations and the zero point of the ephemeris 
(Mor ris"et al. 1990t the cycle count is lost. A clear sinu- 
soidal variation is present; a fit of a sinusoid gives a pe- 
riod of 1.7+0.2 hours, consistent with the period of 1.85 
hours found before .Morris et al. 19901 fMukai et al. 2001. . 
We produced lightcurves in different energy bands to search 
for variations in the amplitude of this modulation (0.3-1.5 
keV, 1.5-10 keV). We fixed the period and the phase of the 
sinusoid to the values obtained from the fit to all energy 
bands combined. The amplitude in the 0.3-1.5 keV band 
was (2.0+0.5) x 10~ 2 counts s _1 whereas that in the 1.5-10 
keV band was (1.3 ± 0.5) xlO -2 counts s _1 (errors are la 
single parameter errors). We note that the 25-30 per cent 
pile-up could in principal have affected the amplitude of the 
sinusoidal modulation as a function of energy. 



3 DISCUSSION 

We observed MS 1603.6+2600 with the ACIS-S instrument 
onboard the Chandra satellite for nearly 7 ksec. Only one 
source was detected and its position is consistent with the 
USNO B1.0 position of the optical counterpart. The source 
was piled-up but still useful spectroscopic parameters could 
be derived by using the pile-up model of |Davis (200Tj\ The 
best-fit spectrum is that of an absorbed power law with 
photon index ~2 and an interstellar absorption consistent 
with the values found by |Mukai et al. (2001)| The power law 
index of 2 is typical for a low-mass X-ray binary (LMXB). 
|Morris et al. (1990)| reported that MS 1603.6+2600 has a 
hard spectrum using Einstein observations. 

The lightcurve shows a sinusoidal modulation with a 
period consistent with that derived previously from opti- 
cal photometric and X-ray observations IMorris et al. 19901 
IMukai et al. 2001.. Similar to the 1991 ROSAT PSPC ob- 
servations ( Hakala et al. 1998J , no partial X-ray eclipse is 
present during our 2002 Chandra observation. However, 
eclipses have been observed during a 1997 ASCA observa- 
tion. |Mukai et al. (2001) | found a clear dependence of the 
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Figure 1. Radial profile of the source. The effect the pile-up has on the radial distribution of the point spread function is clearly visible 
as a reduction of the count rate within the inner pixel. 



Table 1. Best fit parameters of the spectra of MS 1603.6+2600. All quoted errors are at the 90% confidence level. In deriving the value 
for the interstellar absorption, the local absorption due to the Chandra optical blocking niters was accounted for. The first line describes 
the best— fit parameters obtained when using XSPEC whereas on the second line the best— fit parameters are given when using ISIS. 



(xlO 21 cm" 2 ) 


Index 10 — 3 photons keV -1 cm -2 s~ 1,b 


Pile— up parameter 

a c 


Reduced 
X 2 /d.o.f. 


1.5±0.3 
1.2±0.3 


2.05+0.15 2 - 2 -oa 
1.9±0.2 2+3 


0.76±0.25 
0.77+0.41 


1.0/173 
1.1/184 


b 


a PL = power law 
Power law normalisation at 1 keV. 







c Parameter has a hard upper limit of 1. 



lightcurve profile on the photon energy. They observed an 
eclipse at X-rays (<2 keV) whereas at energies >2 keV the 
eclipse is much less pronounced; we find no evidence for a 
change in the amplitude of the sinusoidal lightcurve profile 
as a function of energy. 

Perhaps the fact that we do not find evidence for 
such behaviour ties in with the behaviour observed at 
optical wavelengths. If the overall source luminosity is 
higher, the eclipses in optical are less pronounced; at 
the highest optical luminosities they disappear completely 
HMorris et al. 19901 . Indeed, the 0. 1-10 keV unabsorbed flux 
at the epoch of the observations of |Mukai et al. (2001)| (~5— 
8xl0~ 12 erg cm -2 s" 1 ) is a factor of 1.5 lower than the un- 
absorbed 0.1-10 keV flux we measured in our Chandra ob- 
servations (1.2xl0 -11 ergcm _2 s _1 ). 

As mentioned in the introduction it has been considered 



that MS 1603.6+2600 is a Cataclysmi c Variable (CV) and 
not an LMXB. However, as discussed in |Hakala et al. (1998)| 
and |Mukai et al. (2001)| the optical and X-ray properties of 
MS 1603.6+2600 are difficult to reconcile with a CV sce- 
nario. So, if MS 1603.6+2600 is a CV it must be an unusual 
one IMukai et al. 200111 . There are three possible LMXB sce- 
narios for MS 1603.6+2600 presented in the literature; below 
we will discuss all three of them. 

Interpreting the observed X-ray flare as a type I X— 
ray burst, [Mukai et al. (2001)1 note that the source distance 
must be large (~75 kpc), consistent with earlier estimates 
(30-80 kpc; IMorris et al. 199011 . They further argue that 
based on the variability of the lightcurve in optical and X- 
rays the source differs from an ADC source like 2S 1822- 
371, and hence it is a dipper source like 4U 1916-053. 
However, variability in ADC sources is not new. For in- 
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Figure 2. Upper panel: Chandra ACIS— S X-ray spectrum (0.3-8 keV) of the source MS 1603.6+2600. The best-fit power— law model 
modified by the combined effects of pile— up (Davis 2001) and interstellar absorption and absorption due to contamination of the optical 
blocking filters of the ACIS instrument derived by using XSPEC is overplotted. Lower panel: data minus model residuals. Since the 
residuals are divided by the error the Y— axis is dimcnsionless (the size of the error bars is unity). 



stance, the soft X-ray transient 4U 2129+47 is an ADC 
source in outburst while in quiescence full eclipses are ob- 
served IGarcia 1994l|Nowak, Heinz fe Begelman 2002| . Fur- 
thermore, if MS 1603.6+2600 is a dipper at 75 kpc 
the absolute visual magnitude, My, of the source would 
be ~0, much larger than that of the high inclination 
source in a 50 minute orbit 4U 1916-053 (My = 5.3) 
and even larger than that of the eclipsing source in 
a 3.9 hour orbital period EXO 0748-676 (My = 1.4; 
|van Pa radijs fc McClintock 1995} . The error on the esti- 
mates of these absolute magnitudes will be smaller than ~ 1 
magnitude since the distance estimates for these systems 
were obtained using radius expansion bursts (4U 1916-053, 
ISmale et al. 19881 EXO 0748-676, IGottwald et al. 1986H . 
Such estimates are typically accurate to approximately 
20 per cent tRuulkcrs et al. 2003). Furthermore, the red- 
dening towards these sources is low and well estab- 
lished (4U 1916-053, ISmale et al. 19881 EXO 0748-676, 
ISchoembs fc Zo eschingcr 1990J. Since the absolute optical 
magnitude scales with the binary orbital period and the X- 
ray luminosity {van Paradijs fc McClintock 1995| l, it is un- 
likely that the absolute optical magnitude of the 1.85 hour 
system MS 1603.6+2600 is larger than that of EXO 0748- 
676. On the basis of this, we conclude that it is unlikely that 
MS 1603.6+2600 is a dipper at a distance as large as 75 kpc. 
The X-ray and optical properties of MS 1603.6+2600 are 
consistent with a dipper scenario in which the distance to 
the source is ~15-20 kpc if either the observed X-ray flare 
was not a type I X-ray burst or the luminosity of the burst 
was unusually low. 



|Hakala et al. ( 1998) | favour a less distant (d~0.25-2.7 
kpc) source assuming MS 1603.6+2600 is a soft X-ray tran- 
sient accreting at a low rate in quiescence. Even though the 
X-ray spectrum of several neutron star soft X-ray transients 
in quiescence is well fit by a black body, often a power-law 
component is present in the spectrum as well. The quies- 
cent X-ray spectrum of SAX J1808. 4-3658 was well fit by 
a single power law fCampan a et al. 2002) similar to the X- 
ray spectrum of MS 1603.6+2600 albeit somewhat harder 
(the power-law index for SAX J1808.4-3658 was 1.5). The 
quiescent spectrum of black hole candidate SXTs in qui- 
escence is well fit by a power law with an index of ~2 
{Kong et al. 2002) . 

Several systems containing a neutron star compact 
object accreting at a low rate have been identified 
JCornelisse et al. 200211 . An outburst has never been ob- 
served from these systems. If MS 1603.6+2600 is indeed 
at a distance of approximately 1 kpc then the distance 
modulus is ~10. Given the observed V band magnitude 
of MS 1603.6+2600 (R=19.4, V~19.7; IMorris et al. 19901 
Mv is ~10. Given the fact that many disc lines and 
even the Bowen blend are observed in the spectrum 
JMorris et al. 19901 the disc absolute magnitude must be 
close to this whereas that of the companion star must be 
less than 10; constraining the spectral type of the compan- 
ion star to late type M dwarfs. Such a star indeed fits the 
Roche lobe assuming a neutron star mass of 1.4 Mq. The 
luminosity of MS 1603.6+2600 is close to 10 33 erg s" 1 if the 
distance is 1 kpc. This luminosity is typical for a neutron 
star LMXB in quiescence. Using the same line of reason- 
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Figure 3. Lightcurve of MS 1603.6+2600 in the 0.3-10 keV energy band, excluding data from a circular region with a 2 pixel radius 
centred on the position of MS 1603.6+2600. Each bin is an average of 400 seconds of data. Clearly visible is the sinusoidal variation at 
the orbital period of ~1.85 hours (~6600 seconds). The solid line is the best-fit sinusoid to the data, with a period of 1.7+0.2 hours. 
Time zero is MJD 52344.346. Note that the background is not subtracted. 



ing one can show that a late type M star also nearly fits 
the Roche lobe of a nearby black hole candidate (with an 
assumed black hole mass of 6 Mq) SXT in quiescence. 

However, strong He II and He I lines are present in 
the spectrum of MS 1603.6+2600. Such He features are 
present in the spectra of actively accreting low-mass X- 
ray binaries l |van Paradijs fc McClintock 1995| but not in 
the spectrum of the quiescent soft X-ray transient Cen X- 
4 jvan Paradijs et al. 19871 IChevalier et al. 19891 . In this 
"quiescent soft X-ray transient" scenario the X-ray flare ob- 
served by |Mukai et al. (2001)| cannot have been a type I X- 
ray burst for its luminosity was much too low for a distance 
of ~1 kpc. X-ray flaring in quiescence is not unfeasible since 
optical flares have been found in both neutron star and black 
hole soft X-ray transients in quiescence ( |Hynes et al. 2002] 
|Zurita, Casares fc S hahbaz 2003|> . 

In ADC sources only a fraction of the true source 
luminosity is observed due to scattering in the ADC. In 
case of 4U 2129+47 a type I X-ray burst was observed 
for which the luminosity was a factor of ~500 lower than 
the E ddington luminosity jCarcia fc Grindlay 1987^ . Simi- 
larly, White & Holt (1982)| argue that the intrinsic source 
luminosity of ADC sources should be > 10 37 erg cm -2 s ; 
this is higher than the luminosities one would derive for 
ADC sources if one assumed that the sources are viewed di- 
rectly. Finally, |Jonker fc van der Klis (2001)| presented evi- 
dence based on the observed spin-up of the pulsar in the 
ADC source 2S 1822-371 showing that the intrinsic source 



luminosity is a factor of ~100 higher than the observed lumi- 
nosity. If we assume that we only observe between 1 / 500 and 
1/100 of the intrinsic source luminosity of MS 1603.6+2600 
the distance would be ~ 11-24 kpc assuming that the in- 
trisic burst flux (i.e. 100-500 times the observed flux) of the 
burst in MS 1603.6+2600 observed by ASCA corresponds to 
the Eddington luminosity (we take for the Eddington lumi- 
nosity ~ 2 x 10 38 ergs _1 ; the observed burst peak flux was 
- 3 x 10~ n erg cm" 2 s" 1 : IMiikai et al. 200l")l . Such a dis- 
tance would imply an My intermediate between that of the 
high inclination 50 minute binary 4U 1916-053 and that of 
the 3.9 hour binary EXO 0748-676, as expected. For a dis- 
tance of ~ll-24 kpc the system would be ~8— 17 kpc above 
the Galactic plane. As mentioned before |Mukai et al. (2001)] 
dismiss the identification of MS 1603.6+2600 as an ADC 
source on the ground that MS 1603.6+2600 is different 
from 2S 1822-371 in many facets. However, it is unclear 
if this is enough to dismiss MS 1603.6+2600 as an ADC 
source since the group of ADC sources is highly disparate. 
For instance, 2S 1822-371 harbors a high magnetic field 
neutron star plonker fc van der Klis 2001|l . 4U 2129+47 
is a transient system IGarcia "19941 . and 2S 0921-630 
has a sub-giant companion star (P or b = 9.02 days; 
|Branduardi-Raymont et al. 1983^ . 

We conclude that the classification of MS 1603.6+2600 
depends strongly on whether the X-ray flare observed by 
|Mukai et al. (2001)| was a type I X-ray burst or not. If it 
was MS 1603.6+2600 is most likely an ADC source. If the 



6 P. G. Jonker et al. 



flare was just an X-ray flare and not a type I X-ray burst 
MS 1603.6+2600 is either a nearby soft X-ray transient 
in quiescence or a high inclination LMXB (possibly ADC) 
source located well above the Galactic plane. 
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